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1 Introduction 

The nuclear rotation is an example of a collective motion for which time- 
dependent linear combinations of stationary states of a given spin can be 
identified with a rotating wave packet. A description of such states can 
be performed witin a mean-field theory in which properties of the system 
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are determined by the one-body density matrix. Because the time-reversal 
symmetry is broken for the rotating states, the density matrix, as well as 
the resulting mean field have both time-even and time-odd components. 

Properties of nuclear time-even mean fields are known rather well, be- 
cause they are reflected in multiple static phenomena which can be studied 
experimentally, cf. review On the other hand, very little is known about 
properties of the time-odd mean fields. The present study is devoted to 
an attempt to analyse properties of these time-odd components in rapidly 
rotating superdeformed nuclei. In particular, we aim at studying the phe- 
nomenon of identical bands, which may provide a unique information on 
poorly known sector of the nuclear effective interaction. 

In Sect. |2| we present a construction of the energy-density functional 
based on a few simple assumptions concerning its structure. Section |3| 
presents results of selfconsistent cranking method applied to nuclei in the 
74~150 region; it supplements the results of similar calculations presented in 
Ref. §. 

2 Local density approximation 

The density of nuclear matter in the interior of atomic nuclei has a well de- 
fined value, called the saturation density, which is independent of the nuclear 
size. This property of matter composed of strongly interacting nucleons is a 
basic feature of the nucleon-nucleon interaction which (depending also on the 
intrinsic symmetries of the nucleonic states) gives the strongest attraction 
at a certain distance and becomes weaker or repulsive at other distances. 
Therefore, as a first approximation one can consider a state of the matter 
in a given point of the nucleus as given by that of the infinite matter at the 
saturation density. In other words, properties of nuclear matter in a given 
point weakly depend on other points of the nucleus. 

Such an observation underlies the local density approximation (LDA) 
in nuclear physics It has been adopted from similar ideas developed 
in atomic physics [Q] and reflects fundamental properties of many-fermion 
systems. In such systems, the Pauli correlations play a very important role 
in defining properties of states near its ground state. To some extent, de- 
tails of interactions are less important, and in particular, the local density 
approximation can be valid irrespectively of the interaction range. 

Because of the finite nuclear size, the simplest version of the LDA has to 
be corrected by taking into account not only the density itself but also its lo- 
cal derivatives. This amounts to extending the Thomas- Fermi approximation 
beyond its extreme version, or to applying the systematic Wigner-Kirkwood 



2 



h expansion In a phenomenological approach one may use the ideas of 
the LDA to construct the energy-density functionals, and in what follows we 
proceed along these lines. 

Principal assumptions of such a LDA can be formulated as follows: 

• The total energy of the nucleus is given by the integral of the local 
energy density 7^(r'), 

£ = J d^rn{r), (2.1) 

• The energy density depends on the nuclear one-body density matrix 
and its derivatives up to the second order. 

Neglecting for simplicity the isospin degree of freedom, i.e., considering 
one type of nucleons only, the density matrix in spatial coordinates can be 
defined by 

p{ra,r'a') = (^>|a+(rVXro-)|$), (2.2) 

where r is the position and a is the spin of a nucleon, while |<1>) is a many- 
body nuclear wave function. 

The spin degrees of freedom can be separated by defining the scalar and 
the vector parts of the density matrix, p{r, r') and s(r, r'), respectively, i.e., 



p{r,r') = '^p{ra,r'a) (2.3a) 
s{r,r') = ^p{ra,r'a'){a'\cr\a). (2.3b) 

o-cr' 

Since the density matrix is hermitian, its scalar and vector parts have the 
following properties with respect to exchanging spatial arguments: 

p*{r,r') = p{r',r) (2.4a) 
s*{r,r') = s(r',r). (2.4b) 

On the other hand, the density matrices corresponding to time-reversed 
states read: 

/(r,r') = p*{r,r') (2.5a) 

s^(r,r') = -s*(r,r'). (2.5b) 

These properties mean that for the time-even state the scalar density matrix 
is a real and symmetric function of spatial arguments, whereas the vector 
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Table I: Local densities of a fermion system up to the second order in deriva- 
tives with respect to the relative coordinate 



Type of density 


Order 


Defintion 


Matter: 





p{r) = p{r,r) 


Current: 


1 


j(r) = (l/2i)[(V-V')p(r,r')W 


Kinetic: 


2 


r(r) = [V ■V'p{r,r%=r' 


Spin: 





s{r) = s{r, r) 


Spin current: 


1 




Spin kinetic: 


2 


T{r) = [V- V's(r,0].=.' 



density matrix is imaginary and antisymmetric. Therefore, the local parts 
of the scalar and of the vector densities, 

p{r) = p{r,r) (2.6a) 
s(r) = s(r,r), (2.6b) 

respectively, do not and do change the sign under the time reversal. Hence 
the matter density p{r) is time-even, and the spin density s{r) is time-odd. 

The local derivatives of the density matrix have been defined in Ref. |5| 
and are summarized in Table |. (See Ref. Q for another possible set of 
definitions.) A systematic construction |5| of these derivatives up to a given 
order consists in acting on the scalar and vector density matrices ( ^.3]) with 
the operators (V— V')/2i and setting the arguments equal, r=r', after the 
differentiation is completed. This amounts to calculating derivatives with 
respect to the relative coordinate r—r' at the relative coordinate equal to 
zero. At this point, the derivatives (V+V')/2 with respect to the total 
coordinate r+r' are not considered, because they amount to calculating the 
derivatives after setting r=r', see below. 

In the first order one obtains two new local densities, the current den- 
sity j and the spin current density Jf^uir), Table |. In the second order 
one has to act on the scalar and vector densities with the tensor operator 
(V^— V^)(Vj/— V[,). However, the resulting traceless-symmetric-tensor den- 
sity is not interesting, because there is no other zero-order tensor density 
with which it could have been contracted to construct the corresponding 
term in the energy density, see below. Therefore, only the scalar part of the 
second-order operator should be considered. Due to the identity 

(V-V')^ = (V + V')^-4V- V, (2.7) 



4 



the action of (V — V')^ can be expressed by that of V ■ V', and therefore, 
following Ref. one uses kinetic densities shown in Table |[ 

We may now proceed by constructing the local energy density as a sum of 
terms depending on the above local densities. Such a construction is based 
on the following rules: 

• Every term should be quadratic in local densities. 

An a priori arbitrary coupling constant can be used in front of every 
term. It is well known that the energy densities which are quadratic 
in the matter density do not lead to saturating systems, and therefore, 
a density dependence of the coupling constants should be allowed. 

• The energy density should be invariant w^ith respect to parity, 
time reversal, and rotations. 

This invariance should be considered in a limit of unbroken symme- 
tries. In principle, an extension of the functional beyond this limit is 
arbitrary. For example, for broken parity symmetry one could use in 
the functional either the term C^p^, or two independent terms C^^ p\ 
and C p"^, where p±={p i: p^)/2, and where p^ corresponds to the 
parity inversed state. It is obvious that for the unbroken symmetry 
both choices give the same result, provided that C=C^. Moreover, 
for both choices the functional is invariant with respect to parity. On 
the other hand, for the broken symmetry the first choice gives the 
following contribution to the total energy: 

J d?rCPp^ = J d^rCP{pl + pi), (2.8) 

and therefore it requires that =0^^ , while in the second choice 
these two coupling constants can be chosen independently. These am- 
biguities are even larger for density-dependent coupling constants. Up 
two now, mostly the first choice has been considered, i.e., the function- 
als are constructed in the unbroken symmetry limit. 

• Terms beyond the second order are not taken into account. 

The order of the given term is defined as a sum of orders of both 
densities on which this term depends. The requirement of neglecting 
the terms beyond the second order is consistent with not considering 
the densities beyond the second order, see Table |. This amounts to 
an expansion of the energy density up to the second order in relative 
coordinates. 
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Table II: Construction of all possible terms up to second order in the energy 
density of a fermion system (for comments see text). 



1 


2 


3 


4 


5 


Density 


Derivative 


Order 


Symmetry 


Energy density 








T 


P 


space 


prime 


dependent 


p{r) 







+ 


+ 


scalar 








Vp 


1 


+ 


— 


vector 




(Vp)2 




Ap 


2 


+ 


+ 


scalar 


pAp 






(VxV)(2)p 


2 


+ 


+ 


tensor 






r(r) 




2 


+ 


+ 


scalar 


pr 




jW(r) 


VJ(0) 


1 

2 


+ 
+ 


— 
+ 


scalar 
vector 


(J(0))2 




J(r) 




1 


+ 




vector 




Vp-J 




V J 


2 


+ 


+ 


scalar 


pV ■ J 






V X J 


2 


+ 


+ 


vector 








(VXJ)(2) 


2 


+ 


+ 


tensor 










1 


+ 




tensor 






s{r) 









+ 


vector 






V • s 


1 






scalar 


(V . sf 






V X s 


1 






vector 




(V X S)2 




(VXS)(2) 


1 






tensor 




((VXS)(2))2 




As 


2 




+ 


vector 


s ■ As 




jir) 




1 






vector 


f 


J • V X s 




Vj 


2 




+ 


scalar 








V X J 


2 




+ 


vector 


s -V X j 






(Vxj)(2) 


2 




+ 


tensor 






T{r) 




2 




+ 


vector 


s-T 





The construction of the energy density is illustrated in Table |Tl[ The first 
column gives the densities presented in Table |. The spin current density 
J^u is split into the standard scalar, vector, and tensor parts (J^°\ J, and 
J^^-' , respectively) according to 
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where 

(2.10a) 
(2.10b) 

Y.Jk.- (2.10c) 

The second column of Table ^ presents derivatives of the densities from 
the first column. The derivatives are constructed in the vector-coupled form. 
For example, from the (vector) gradient operator V and the scalar density p 
one can obtain the vector density Vp, the scalar density A/?, and the tensor 
density (VxV)^^V- Since p itself is a zero-order density, the density Vp 
is of the first order and the other two densities are of the second order, as 
indicated in the third column of the Table. 

The derivatives are calculated only up to the second order. Moreover, 
the rank-3 tensors are not taken into account and the derivatives of J^J (r) 
are not included in the Table, because there are no other densities with 
the same symmetries to construct terms in the energy density. Apart from 
these two omissions, the Table presents a complete set of densities and their 
derivatives which can be obtained within the specified limits. 

The fourth column of the Table gives the symmetries of every density in 
the unbroken symmetry limit, i.e., assuming that p is a time-even, parity- 
even scalar density, and s is a time-odd, parity-even vector density. Then 
the symmetries presented in the Table result from the fact that the gradient 
operator V is a time-even, parity-odd vector operator. 

Finally, the fifth column of the Table presents the terms in the energy 
density which can be constructed from the densities listed in the first two 
columns. In each case a term is obtained by forming a scalar, time-even, 
and parity-even product of two densities, i.e., by multiplying densities hav- 
ing exactly the same symmetries. Two avoid the double-counting, a given 
density is multiplied only by densities appearing higher in the Table. The 
rule of not including terms obove the second order is also respected. 

For example, according to this construction scheme, the matter density 
p gives the simplest term in the energy density, p^, as shown in the first line 
of the Table. Similarly, the derivative V/j gives rise to the term (Vp)^. This 
term is listed in the column of "dependent" terms, because for a density- 
independent coupling constant it can be expressed by the "prime" term 
pAp given in the third line. This can be done by integrating by parts the 



jj,!/ 



UflJ 
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integral of the energy density (2.1). The corresponding expressions for five 
"dependent" terms follows: 



{Vpf = -pAp, (2.11a) 

Vp-J = -pV-J, (2.11b) 

{Vxsf = -{V-sf-s-As, (2.11c) 

E((Vx4l))^ = |(V•.)^ (2.11d) 

jf-Vxs = s-Vxj. (2. lie) 



Altogether, this construction procedure gives eighteen terms in the en- 
ergy density, of which thirteen are "prime" and five are "dependent" . Three 
terms correspond to the squares of the spin current density J^u{r) for the 
three different values of the intermediate-coupling angular momentum. Usu- 
aly, these three independent terms are considered together in the following 
single term: 

with a common coupling constant. 

It turns out that each density, or a derivative of the density, gives rise to 
at most one "prime" term in the energy density, and these terms and coupling 
constants can be conveniently labelled by the corresponding densities. We 
then obtain the following expression for the interaction energy density: 

7i^^\r) = n'^^^'^ir) + n'^'^'^ir), (2.13) 

with 

^even^^) = C p^ + C^" pAp + pT + C^^ T 

+C^-^pV ■ J, (2.14a) 
n'"^'^{r) = C"s^ + C^'s- As + C^s-T + C^f 

+C^^s-{V xj) + C^'{V -sf. (2.14b) 

Each of these terms can have the isoscalar and the isovector form, and the 
kinetic energy density has to be added to obtain the total energy density of 
a nuclear system. 

Apart from the last term in 7i°'^'^{r), all these terms appear in the energy 
density corresponding to the Skyrme interaction |2|. The additional term 
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in Eq. (2.14b) can only be obtained if the tensor interaction is added to the 
Skyrme force, i.e., when the results of Ref. Q are generalized to the case of 
broken time-reversal symmetry. In the following we do not consider such a 
possibility, and we set C^*=0. 

By varying the total energy ( p.!]) with respect to the single-particle wave 
functions one obtains the Hartree-Fock time-even and time-odd mean fields 
M. In the convention corresponding to the energy density of Eq. (|2l|), 
the relevant expressions are given in Ref. |Q| and will not be repeated here. 
Instead, in the following section we present the results of the Hartree-Fock 
cranking method applied to the description of superdeformed bands. 

3 Superdeformed bands in ^^^Dy, ^^^Tb, and ^^^Gd 

In Fig. H we show the results of calculations for the yrast superdeformed 
band in ^^-^Dy and for the first excited band in ^^^Tb, denoted by ^^^Dy(l) 
and ^^^Tb(2), respectvely. Details of the calculation methods can be found 
in Ref. [^. The coupling constants of the energy-density functional have 
been determined from the parameters of three Skyrme forces, SkM* P], SkP 



[P, and SIII |10|. The left-hand-side panels show the results obtained for 
the complete functionals, while those at the right-hand-side correspond to 
omitted time-odd terms TC°'^'^{r)=0, i.e., to 

C = C^' = C^ = = C^^ = 0. (3.15) 

For the ^^^Th{l) and i^°Gd(2) bands, the analo gous results are shown in 
Fig. I 

The omission of time-odd terms in the Hartree-Fock approach is analo- 
gous to the standard phenomenological rotating mean-field approach where 
the selfconsistent changes of the mean field due to rotations are not taken 
into account. Results of latter type calculations, for the same SD bands as 



considered in the present study, can be found in Refs. |11| and |12| . 

For the complete Skyrme functionals, the dynamic moments JT^^^ of the 
^^^Dy(l) band obtained for SkM*, SkP, and SIII interactions are very similar, 
see part (a) in the left-hand-side panel of Figs. |l| and |2[ This is not the 
case when the time-odd terms are omitted, as show in the corresponding 
right-hand-side panels. For SkM* and SIII, the omission of the time-odd 
terms significantly decreases the values of J'^'^'^ , while for SkP the changes of 
are much smaller. This feature is related to the values of the isoscalar 
nuclear- matter effective mass, which for SkP is equal to the free nucleon 
mass, m*=m, while for SkM* and SIII is smaller by the factor of 0.79 and 
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0.76, respectively. As a consequnce, for SkP the isoscalar coupling constant 
Cq is equal to zero and the modifications introduced by omitting the 
time-odd terms, Eq. (|3.15| ), do not influence J'^'^^ in a very strong way. On 
the other hand, whenever this coupling constant is relatively large, as is for 
SkM* and SIII, the corresponding changes in are well prononced. 

These results illustrate the importance for the rotational properties of 
the nucleus of the term in the energy density. The cranking term mainly 
induces the nonzero flow of nuclear matter, as given by the current density j, 
whch influences the time-odd nuclear mean field provided the corresponding 
coupling constant is nonzero. 

In spite of the strong influence on the values of J'^'^\ the omission of 
the time-odd terms has a very small effect on the sameness of ^'-^^ in pairs 
of identical bands. In parts (b) of Figs. || and § we show the the relative 
dynamical moments 5J^'^\ i.e., the differences between the values of JT''-^-* 
calculated for two pairs of identical bands ^^^Th{2)-^^^T)y{l) and ^^°Gd(2)- 
^^^Tb(l). One can see that the omission of the time-odd terms influences the 
values of J^'^^ in a very similar way for both members of each pair. Even if the 
values may change by as much as 15 Ti^/MeV, the relative values are always 
well below 2 Tt^/M.eM . This shows that the sameness of J^'^'^ is not governed 
by the time-odd terms in the mean-fields, but rather can be attributed to 
general geometric properties of the uderlying orbital 7r[301]l/2(r=+i). 

A different patterns is obtained for relative alignments shown in parts 
(c) of the figures. Here the omision of the time-odd terms leads to important 
changes of the relative alignments, and also different relative alignments are 
obtained of the three Skyrme forces studied. 

In Figs. |3| and ^ we present a detailed analysis of the time-odd terms ap- 
pearing in the energy density for the SkP interaction. As shown in the 
legends, the results presented in the left-hand-side panels correspond to 
the following conditions: (i) the complete functional, (ii) Cf=Cf=0, (iii) 
C7/^^=C7f =C/=0, and then (iv) Q=C7/^^=Cf =C/=0. In the right-hand- 
panels, in addtion to the last of these conditions we set to zero either one, 
or both of the and C^^ coupling constants. 

As discussed above, the value of the effective mass m*=m renders the 
dynamical moments J'^'^\ the relative dynamical moments, and the relative 
alignments 51 almost independent of the coupling constants. (A weak 
residual dependence results from a non-zero value of the isovector coupling 
constant C;^, see Ref. |^].) On the other hand, the values of the relative 
alingments strongly depend on the C^-' coupling constants, and are closest 
to the experimental data when these coupling constants are equal to zero. 
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4 Conclusions 



In the present study we have presented a systematic construction of the 
energy-density functional in the frame of the nuclear local density approx- 
imation. Such a construction is based on a few simple prescriptions and 
leads to the functional identical to that obtained by using the Skyrme ef- 
fective interaction. The only exception is the term given by the square of 
the divergence of the spin density, which can only be obtained if the tensor 
component is added to the Skyrme interaction. 

Following Ref. we have here presented the analysis of the time-odd 
components of the mean field of rotating superdeformed nuclei. Special 
attention has been devoted to properties of the SkP Skyrme interaction 
which is characterized be the effective mass m*=m. In particular, we have 
shown that the sameness of the dynamical moments J^"^^ does not depend on 
whether the time-odd terms are, or are not taken into account. On the other 
hand, the sameness of the corresponding alignments does depend on these 
time-odd terms. Therefore, one may expect that a systematic investigation of 
the rotational alingments may serve as a tool for establishing the properties 
of time-odd component of the rotating nuclear mean field. 

References 

[1] S. Aberg, H. Flocard, W. Nazarewicz, Annu. Rev. Nucl. Part. Sci. 40, 
439 (1990); and references therein. 

[2] J. Dobaczewski, J. Dudek, Phys. Rev. C52, No. 4 (1995). 

[3] P. Ring, P. Schuck, The Nuclear Many-Body Problem, Springer- Verlag, 
1980. 

[4] R.M Dreizler, E.K.U. Gross, Density Functional Theory, Springer- 
Verlag, 1990. 

[5] Y.M. Engel, D.M. Brink, K. Goeke, S.J. Krieger, D. Vautherin, Nucl. 
Phys. A249, 215 (1975). 

[6] K.-H. Passler, Nucl. Phys. A257, 253 (1976). 

[7] Fl. Stancu, D.M. Brink, H. Flocard, Phys. Lett. 68B, 108 (1977). 

[8] J. Bartel, P. Quentin, M. Brack, C. Guet, H.B. Hakansson, Nucl. Phys. 
A386, 79 (1982). 



11 



[9] J. Dobaczewski, H. Flocard, J. Treiner, Nucl. Phys. A422, 103 (1984). 

[10] M. Beiner, H. Flocard, N. Van Giai, P. Quentin, Nucl. Phys. A238, 29 
(1975). 

[11] W. Nazarewicz, R. Wyss, A. Johnson, Nucl. Phys. A503, 285 (1989). 
[12] I. Ragnarsson, Nucl. Phys. A557, 167c (1993). 

Figure 1: Calculated dynamical moment ^^^^ for the yrast band of ^^^Dy, 
part (a), the relative dynamical moment dj^"^^ calculated for the ^^^Tb(2) 
and ^^■^Dy(l) bands, part (b), and the relative alignment 51 between these 
two bands, part (c). Left-hand-side panels show the results for complete 
Skyrme functionals of the SkM*, SkP, and SIII interactions, and the right- 
hand-side panels show the results with omitted time-odd terms. The exper- 
imental points are denoted by asterisks. Note the scale in (b) expanded five 
times as compared to (a). 



Figure 2: Same as Fig. |l], but for the ^^^Tb(l) band, (a), and for the differ- 
ences between the ^^'^Gd(2) and -'^^-'^Tb(l) bands, (b) and (c). 

Figure 3: Similar as in Fig. |l|, but for different time-odd terms omitted in 
the energy density, see text. 

Figure 4: Same as Fig. ^ but for the ^^^Tb(l) band, (a), and for the differ- 
ences between the ^^°Gd(2) and ^^^Tb(l) bands, (b) and (c). 
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